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Abstract 

We investigate the B and D mesons in the relativistic quark model by applying 
the variational method with the Gaussian wave function. We calculate the Fermi 
momentum parameter p F , and obtain p F = 0.50 ~ 0.54 GeV, which is almost 
independent of the input parameters, a s , m&, m c and m sp . We then calculate the 
ratio Jb/ Id, and obtain the result which is larger, by the factor of about 1.3, 



than tJMd/Mb given by the naive nonrelativistic analogy. This result is in a 
good agreement with the recent Lattice calculations. We also calculate the ratio 
(Mb* — M B )/(M D * — M D ). In these calculations the wave function at origin ^(0) 
is essential. We also determine p F by comparing the theoretical prediction of the 
ACCMM model with the lepton energy spectrum of B — > evX from the recent 
ARGUS analysis, and find that p F = 0.27 ± 027 GeV, when we use m c = 1.5 
GeV. However, this experimentally determined value of p F is strongly dependent 
on the value of input parameter m c . 
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1. Introduction 



The B meson physics is important in the present high energy physics, since 
it gives us the information on the CKM matrix elements, V c b and V u b, and it is 
expected to provide the CP violation phenomena. In order to extract the value of 
V u b from the B meson decay experiments, the method of separating the B — ► X u lv 
events from the B — > XJu ones at the end point region of the lepton energy 
spectrum has been used Jl|, |2[ . On the other hand, the method of using the hadronic 
invariant mass spectrum has also been suggested recently [||. For the analysis of 
the inclusive semileptonic decay process, the ACCMM model [|J has been most 
popularly imployed, where the Fermi momentum parameter p F is introduced as the 
most important parameter. The value of p F ~ 0.3 GeV has been commonly used 
for the experimental analyses without clear experimental or theoretical supports. 
In reference H we calculated p p in the relativistic quark model using the variational 
method with the Gaussian wave function, and obtained the result of p F =0.54 GeV. 
In this paper we also study the ratios, /b//d an d (Mb* — M B )/(M D * — Mo), using 
the same method. 

When one treats the heavy-light meson in analogy with the nonrelativistic 
situation, one expects /b//d — \Jm~oJMb, since the reduced masses of the light 
and heavy quark systems of the B and D mesons have similar values, and fpMp = 
12 1-0(0) | 2 by the van Royen-Weisskopf formula for the pseudoscalar meson P ||, 
where ip(0) is the wave function at the origin of the relative motion of quarks. If one 
uses the relation /b//d — J M D / M B with the supplementary relation f D / f Ds ~ 
\Jmd/m s , one can obtain the value of fs from that of fo s , which is rather well 
established from the branching ratio of the B° — > D + D~ or — > process ]7). 
However, our calculation of fp in the relativistic quark model, which we present 
in this article, shows that the above consideration with nonrelativistic analogy 
is much deviated by the relativistic motion of the light quark in the heavy-light 
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pseudoscalar meson P. This deviation has also been exposed by the recent Lattice 
calculations, since they give rather close values for f B and f D 0, [j| . This situation 
can be understood clearly within our relativistic calculation. 

The potential model has been successful for ip and T families with the nonrel- 
ativistic Hamiltonian, since their heavy quarks can be treated nonrelativistically. 
However, for D or B meson it has been difficult to apply the potential model 
because of the relativistic motion of the light quark in D or B meson. In our cal- 
culation we work with the realistic Hamiltonian which is relativistic for the light 
quark and nonrelativistic for the heavy quark, and adopt the variational method. 
We take the Gaussian function as the trial wave function, and obtain the ground 
state energy (and the wave function) by minimizing the expectation value of the 
Hamiltonian. Using the Gaussian wave function calculated as above, we get the 
wave function at the origin ip(0), and with which we can calculate the decay con- 
stant of the heavy-light pseudoscalar meson from the van Royen-Weisskopf formula. 
Through this procedure we obtain the ratio We also calculate the ratio 

{Mb* — Mb) /(Md» — Md) from the chromomagnetic hyperfine splitting formula, 
where the information of ip(0) is essential. Finally, we compare the value of the 
Fermi momentum p F given in our calculation with the lepton energy spectrum 
data of the semileptonic decay process, and find that it is just outside off one a 
standard deviation. 

In Section 2, using variational method in relativistic quark model we calculate 
the value of the parameter^, and the ratios and (Mb* —Mb)/ {Mb* — Mb)- 

We also extract p F by comparing the theoretical prediction of the ACCMM model 
with the whole region of electron energy spectrum of B — > evX in Section 3. Sec- 
tion 4 contains the conclusions. 
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2. Variational method in relativistic quark model, and calculations 
of f B /f D and (M B * - M B )/(M D * - M D ) 



For the B meson system we treat the 6-quark nonrelativistically, but the u- or 
ci-quark relativistically with the Hamiltonian 

P 2 



H = M + 2M + V p2 + m2 + V(r) ' (1) 

where M = or m c is the heavy quark mass and m = m sp is the u- or d- 
quark mass, i.e. the spectator quark mass in the ACCMM model. We apply the 
variational method to the Hamiltonian ([!]) with the trial wave function 

^(r) = (^)3/V^ 2 / 2 , (2) 
V 71 " 

where \i is the variational parameter. The Fourier transform of ip(r) gives the 
momentum space wave function x(p)> which is also Gaussian, 

x(p) - yhfr^- (3) 

The ground state is given by minimizing the expectation value of H, 

(H) = ty\H\ij>) = E(ji), ^-E{p) = at f, = ft, (4) 

dfi 

and then E = E{jj) approximates tub, and we get jx = p F , the Fermi momentum 
parameter in the ACCMM model. The value of \i or p F corresponds to the measure 
of the radius of the two body bound state, as can be seen from the relation, 

( r \ = 2 1 Qr / 2 U = 3 1 
N ' V 71 " M 2 fj, 

In Eq. (HD, we take the Cornell potential which is composed of the Coulomb 
and linear potentials, 

OL 

V(r) = + Kr. (5) 

r 

For the values of the parameters a c (= |a s ), the quark masses m b and m c , and 
K, we use the following two sets of parameters. The set {A} is the Hagiwara et 



a/.'s [ID | which has been determined by the best fit of (cc) and (bb) bound state 



spectra. And the set {B} is chosen to have the running coupling constants at 
the mass scales of mg and mjj, and the quark masses m& and m c which were 
determined to give the best ip and T masses for the variational ground states. 

{A} a c = 0.47, m b = 4.75 GeV, m c = 1.32 GeV, K = 0.19 GeV 2 , (6) 
{B}af = 0.32, m b = 4.64 GeV; af = 0.48, m c = 1.33 GeV; K = 0.19 GeV 2 . (7) 

With the Gaussian trial wave function (|j) or @, the expectation value of each 
term of the Hamiltonian (p]) is given as follows: 

(^7> = (X(P)I7777IX(P)> = -777V 2 



2M' ^"2^'^" 4M^ ' 

1// /' x 

7T JO 



p 2 + m 2 ) = (x(p)I\/p 2 + m2 |x(p)) = "~ 7= / e x2 \l x<2 + {m/n) 2 x 2 dx 



rv 2 
(y(r)) = (^(r)|-^ + Kr |^(r)) = — (- acfi + K/fj). (8) 

r ^/7^ 



Then we have 



E(ji) = (H) (9) 
= M + -77/i 2 + ^=(— a c /i + K//j) + ^= ( e' x \/x 2 + (m//i) 2 x 2 c/x. 

4M V 7T V7T J0 



For more details on this procedure of the variational method, see Ref. ||. 

With the input value of m = m sp = 0.15 GeV, which is the value commonly 
used in experimental analyses, we minimize E(/a) of @, and then we obtain for 
the B meson, 

p F (B)=fL = 0.54 GeV, E(B) = 5.54 GeV for {A}, (10) 
p, = 0.50 GeV, E(B) = 5.52 GeV for {B}. 

The B meson mass is lowered from the above values if we include chromomagnetic 
hyperfine splitting corrections. For comparison, let us check how much sensitive 
our calculation of p F is by considering the case where m = m sp = 0. For m sp = 
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the integral in (|9|) is done easily and we obtain the following values of // = p F , 

ji = 0.53 GeV, E(B) = 5.52 GeV for {A}, (11) 
Ji = 0.48 GeV, E(B) = 5.49 GeV for {B}. 



As we see in Eq. (|11]), the results are similar to those in (|I0|), where 



0.15 GeV. We expected this insensitivity of the value of p F on m sp because the 
value of m sp , which should be small in any case, can not affect the integral in 
(H) significantly. We also note that the theoretically determined value of p F (B) 
is completely independent of the input value of m c , as can be seen from Eq. (||). 
Following the same procedure, we next obtain the results for the D meson with 
m sp = 0.15 GeV, 

p F (D)=fi = 0.45 GeV, E(D) = 2.21 GeV for {A}, (12) 
Pl = 0.46 GeV, E(D) = 2.21 GeV for {B}. 

The decay constant fp of a pseudoscalar meson P is defined by the matrix 
element (0|A M |P(g)): 

(0\A lt \P{q)) = iq li f P . (13) 

By considering the low energy limit of the heavy meson annihilation, we have the 
relation between fp and the ground state wave function at the origin ipp(0) from 
the van Royen-Weisskopf formula including the color factor |], [/J : 

fp = ^-\M0)\\ (14) 



where Mp is the heavy meson mass. From (14) we have the ratio of fp and fp,: 



f D y M B 

For the Gaussian wave function of Eq. @, we have 



fn jMn x hM0)l ^ (15) 



Mo) = (^S) 3/2 , (16) 

'7T 



then using the values of p F in ([[(]) and (|12"1), we obtain 



|Wt X 3/2 ^°- 59Xl - 3 -°- 77 for {A}, ( 17 ) 

= 0.59 x 1.13 = 0.67 for {B}. 



From (|T7|) we see that f B j ' f D is enhanced, compared with y Mo /Mb, by the factor 
1.31 for the parameter set {A}, and by 1.13 for the set {B}, which are given by 
the factor of \ip B (0)/ip B (0)\- Sometimes this factor has been approximated to be 



1, and the relation f B j f D ~ ^M B /M B has been used, by treating it in analogy 
with the nonrelativistic case [|?J. However, our calculation shows that this factor 
is indeed important and different from 1 significantly. The factor 1.31 obtained in 
( |T7D for the parameter set {A} is in a pretty good agreement with 1.40 of Ref. M 
and 1.39 of Ref. || of the recent Lattice calculations. 

The mass difference between the vector meson P* and the pseudoscalar meson 
P is given rise to by the chromomagnetic hyperfine splitting: 

V hf = -^s 1 .s 2 X7\-^), (18) 
3m Q m q r 

where rriQ is the heavy quark mass as we used before (i.e. m b or m c ), and m q is 

the constituent quark mass of the light quark, which is the effective mass for the 

baryon magnetic moments [fTI ], 



m u = m d = 0.33 GeV, m s = 0.53 GeV. (19) 
Then the mass difference between B* and B mesons is given by 

Using the values of p F in flT0| ) and (|T2D, we obtain the ratio of M B * — M B and 

M fl » - M D as 

M B , - M B = mc / |^b(0)| \ 2 = rric / p F (B) y 

= 0.28 x 1.73 = 0.48 for {A}, (21) 
= 0.29 x 1.28 = 0.37 for {B}. 



The experimental value [0 is about 0.33, which is larger than the nonrelativistic 
value 0.28 or 0.29, but smaller than our calculated value 0.48 or 0.37. Our cal- 
culated results are at least not much worse than the nonrelativistic values, even 
though there exists somewhat large discrepancy compared to experimental result. 
And this suggests that the behind reason for this discrepancy could be more subtle 
than the nonrelativistic consideration. 

3. Determination of p F from the experimental spectrum 

Until now, we have discussed the theoretical determination of p F in the rel- 
ativistic quark model using the variational method, and its implications to the 
heavy meson masses and the decay constants of the heavy mesons. Now we would 
like to determine the Fermi momentum parameter p F by comparing the theoretical 
prediction with the experimental charged lepton energy spectrum in semileptonic 
decays of B meson. 

As discussed, the simplest model for the semileptonic S-decay is the spectator 
model which considers the decaying 6-quark in the B meson as a free particle. The 
spectator model is usually used with the inclusion of perturbative QCD radiative 
corrections [ 13| . Then the decay width of the process B — > X q lv is given by 



r B(B x qlv) * r,(6 = i^tf^/Qn - (22) 

where m q is the mass of the final g-quark decayed from 6-quark. As can be seen, 
the decay width of the spectator model depends on mjj, therefore small difference 
of nib would change the decay width significantly. 

Altarelli et al. |4[] proposed for the inclusive B meson semileptonic decays their 
ACCMM model, which incorporates the bound state effect by treating the 6-quark 
as a virtual state particle, thus giving momentum dependence to the 6-quark mass. 
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The virtual state 6-quark mass W is given by 

W 2 (p) = ml + m% - 2m By J^ + m% (23) 

in the B meson rest frame, where mg is the B meson mass, m sp is the spectator 
quark mass, and p is the momentum of the spectator quark inside B meson. 

For the momentum distribution of the virtual 6-quark, Altarelli et al. considered 
the Fermi motion inside the B meson with the Gaussian momentum distribution, 

0(p;p F ) = 4tt| X (p)| 2 = -i^e-^, (24) 

where the Gaussian width, p F , is treated as a free parameter. Then the lepton 
energy spectrum of the B meson decay is given by 

(IT B rVmax (IT), 

-TEr{p F ,™> ap ,m q ,m B ) = dp p (p(p;p F ) -jw{m b = W,m q ,m sp ), (25) 

OlEji J CLJi'l 

where p ma x is the maximum kinematically allowed value of p = |p|. The ACCMM 
model, therefore, introduces a new parameter p F for the Gaussian momentum 
distribution of the 6-quark inside B meson, instead of the 6-quark mass of the 
spectator model. In this way the ACCMM model incorporates the bound state 
effects and reduces the strong dependence on the 6-quark mass in the decay width 
of the spectator model. The Fermi momentum parameter p F is the most essential 
parameter of the ACCMM model, as we see in the above. However, the experi- 
mental determination of its value from the lepton energy spectrum has been very 
ambiguous, because various parameters of ACCMM model, such as p F , m q and 
fitted all together from the limited region of end-point lepton energy 
spectrum, and because the perturbative QCD corrections are very sensitive in the 
end-point region of the spectrum. 



Recently, ARGUS [[LJ] extracted the model independent lepton energy spec- 
trum of B — > X c lu for the whole region of electron energy, but with much larger 
uncertainties, as shown in Fig. 1. We now compare the whole region of experimen- 
tal electron energy spectrum with the theoretical prediction of ACCMM model, Eq. 
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fl25|), using p F as a free parameter. We fixed m sp = 0.15 GeV and m q = m c = 1.5 
GeV, which are the values commonly used by experimental analyses. We derive 
the value of p F using \ 2 analysis, and we obtain 

p F = 0.27 ± GeV. (26) 

The minimum y 2 equals to 0.59 with p F = 0.27 GeV. However, the result, Eq. 
(|26D, is found to be strongly dependent on the input value of m c : if we instead 
use smaller m c , both the best fit value of p F and the minimum \ 2 increase, and 
vise versa. In Fig. 1, we also show the theoretical ACCMM model spectrums with 
p F = 0, 0.27, 0.49 GeV, corresponding to dashed-, full-, dotted-line, respectively. 
The experimental data and the theoretical predictions are all normalized to the 
semileptonic branching ratio, BR(B — > evX) = 9.6%, following the result of 
ARGUS @. 

In Section 2, we theoretically calculated the Fermi momentum parameter p F , 
and obtained p F = 0.50 ~ 0.54 GeV. We note that the theoretically calculated 
values are slightly outside off one a standard deviation comparing with the best fit 
value of the experimental data. However, since the experimental spectrum has still 
large uncertainties, we cannot exclude the validity of relativistic quark model for 
calculation of p F yet, nor ACCMM model itself to apply to experimental analyses 
to find CKM parameters, |V^,| and/or |V^/V^|. In near future, once we get much 
more data from asymmetric B factories, it would be very interesting to extract the 
precise value of p F once again. 

In our previous work [|], we investigated the dependence of |T4&/K*| on p F in 
ACCMM model, and found rather strong dependence as a function of a parameter 
P F , 

10 2 x \V ub /V cb \ 2 = 0.57 ±0.11 (ACCMM with p F = 0.3 [15]), 
= 1.03 ±0.11 (ACCMM with p F = 0.5), 
= 1.02 ±0.20 (Isgur et al. [0). (27) 
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As can be seen easily, those values between ACCMM model with p F = 0.3 GeV 
and Isgur et al. model are in large disagreement. However, if we use p F = 0.5 
GeV, the result of the ACCMM model becomes 1.03, and these two models are in 
a good agreement for the value of | V^&/V^b| . 

4. Conclusions 

We conclude that the value p F ~ 0.3 GeV, which has been commonly used in 
experimental analyses, has no theoretical or experimental clear justification, even 
though there has been recently an assertion that the prediction of heavy quark 
effective theory approach |17]], far from the end-point region, gives approximately 
equal shape to the ACCMM model with p F ~ 0.3 GeV. Therefore, it is strongly rec- 
ommended to determine the value of p F more reliably and independently, when we 
think of the importance of its role in experimental analyses. It is particularly im- 
portant in the determination of the value of | / | . A better determination of p F 
is also interesting theoretically since it has its own physical correspondence related 
to the Fermi motion inside B meson. In this context we calculated theoretically the 
value of p F in the relativistic quark model using quantum mechanical variational 
method. It turns out that p F = 0.50 ~ 0.54 GeV, which is not far from the value of 
p F determined by comparing ACCMM model prediction and the model indepen- 
dent lepton energy spectrum of ARGUS measurement, p F = 0.27±{];27 GeV. The 
theoretically determined value of p F is almost independent of input parameters, 
a s , rat,, Tn c , m sp and etc. On the other hand, the experimentally determined value 
is strongly affected by the value of the input parameter m c . 

By using the same framework, we then calculated the ratio /b//d; and obtaind 
the result which is larger, by the factor of about 1.3, than J M^/Mb given by the 
naive nonrelativistic analogy. This result is in a pretty good agreement with the 
recent Lattice calculations. We also calculated the ratio {Mb* —M b ) / (M D * —M D ), 
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whose results suggest the subtlety of the mechanism which gives rise to the value 
of this ratio. 
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Figure Captions 



Fig. 1 The normalized lepton energy spectrum of B — > XJu for the whole region 
of electron energy from the recent ARGUS []I4J measurement. Also shown are the 
theoretical ACCMM model predictions, Eq. (0), using p F = 0, 0.27, 0.49 GeV, 
corresponding to dashed-, full-, dotted-line, respectively. The minimum \ 2 equals 
to 0.59 with p F = 0.27 GeV. We fixed m sp = 0.15 GeV and m q = m c = 1.5 GeV. 
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